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Abstract 

Using a sample of 38.5 million underground muons collected by the MACRO detector we have performed 
an all-sky search for pointlike sources producing excesses of muons above the expected background. The 
d.c. muon flux upper limit at the Earth coming from selected sources is of the order of 10~ 13 cm~ 2 s _1 or 
less. Futhermore we discuss searches for possible modulated and variable sources of muons using different 
j techniques. 
ON \ 

1 Introduction: 

In this work we present the results on the study of muon arrival directions, as seen by the MACRO detector 
(Ahlen et al., 1993), looking for excesses above the expected background in every sky direction. We have also 
searched for steady and pulsed signals from potential cosmic ray sources. 

Features of the MACRO detector are an excellent angular resolution and a large collecting area, when 
compared to other deep underground detectors. Because of these properties, MACRO can be used as a good 
muon telescope. The interest in this kind of analysis arises from the old observation of underground muon 
excesses from the binary system Cygnus X-3, as reported from the SOUDAN and NUSEX collaborations 
_ \ (Marshak et al., 1985; Battistoni et al., 1985). On 1991 January, SOUDAN2 reported the observation of 
unmodulated muon signals from Cyg X-3, correlated with the radio-burst occurring in the same period. These 
^ ■ muon excesses had fluxes of 7.5 x 10~ 10 cm~ 2 s^ 1 (Marshak, 1993). Recently, no experiment detected signals 
O ! from this kind of source, but several observations of X and 7-bursts were reported from other celestial objects. 

During 1997, several 7-bursts from AGN objects, like Markarian 421 and Markarian 501, were reported by 
the Whipple and Hegra experiments (Punch et al., 1992; Quinn, 1997). 

x : 

^ ■ 2 Data Selection and Reduction: 

Oh. The data set used for this analysis is from a long period of acquisition, from the first MACRO run (February 

1989) until January 1999. In order to optimize the track reconstruction and to eliminate possible detector 
malfunctions, strict selection criteria were defined. After all cuts, we selected a sample of 38.5 million events 
I • »-j . (single and doubles), collected over 55248 h of live time. 

2.1 Angular Resolution and Background Simulation: The MACRO pointing precision was cal- 
culated from the distribution of the double-/i angular deviation. We assume that these events at the Earth's 
surface are parallel muon pairs. Then, the overburden rock produces their deviation at the MACRO level 
through multiple Coulomb scattering processes. The space angle 9, containing 68% of the events, defines the 
detector angular resolution. We found 8 = 0.8°, consistent with the value derived from the Moon shadowing 
of primary cosmic rays (Giglietto, 1999). 

For the present search, we assume that our signals are muons produced in atmospheric photoproduction 
processes and the background (or noise) is produced by primary interactions with the terrestrial atmosphere. 
The background has been determinated by a Monte Carlo simulation, assuming an isotropically distributed 
cosmic ray flux. For each run we generated the background events by coupling the arrival direction of each 
muon with 25 arrival times, randomly extracted from the same run. 

3 All-sky d.c. Survey: 

In the search for point-like sources we first started an all-sky survey, without a priori assumptions about 
the source locations; then, we examined selected celestial objects known as potential cosmic ray sources. 



Therefore, we examine sky regions looking for deviations from the expected background larger than random 
fluctuations. 

3.1 Search for d.c. Muon Excesses: As a first step, we generated the sky map with muons collected 
by MACRO by dividing the sky into bins of equal solid angle (AS1= 2. 1 x 10~ 3 sr; Aa = 3°, Asin<5 = 0.04). 
These bins have the same Afi as a nar- 
row cone of half-angle 1.5°. Then, assum- 
ing that a potential source is located at the 
center of the bin, we obtained a sky sur- 
vey in absolute astronomical coordinates 
a and sim5 (map(l) in Table 1). However, 
if the source is near a bin edge, its sig- 
nal will be spread between adjacent bins. 
Then, to reduce the possibility of missing 
sources close to the bin edges, three other 
surveys were also done: map(2) shifted 
the grid one-half bin in a, map(3) shifted 
the grid one-half bin in sine) and map(4) 

which contains the map(2) and map(3) shift. For each solid angle bin in the four surveys, we calculated the 
deviation from the mean in units of standard deviations: 
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map(3) 


-0.22 xl0~ 3 


1.01 
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4(+) 4(-) 


map(4) 


0.3xl0" 2 


1.02 


34/35 


5(+) 4(-) 



Table 1 : Best-fit parameters of the Gaussian fit for the four sky 
surveys. The last columns gives the number of positive and nega- 
tive deviations larger than 3.2a. 
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where N b s (i) is the observed number of events in the bin and N exp (i) is the number of events expected from 
the background simulation. Bins having less than 20 events were removed. 

If the observed muons are distributed according to Gaussian statistics, the a(i) distributions will be a 
Gaussian curve with zero mean and unit standard deviation. Table 1 shows the Gaussian best fit for the four 
surveys: the means are close to zero and the half widths near one. The best-fit parameters and the number of 
fluctuations larger than 3.2<r are also given. We found the larger positive and negative deviations are in similar 
number. We conclude that all deviations are consistent with random background fluctuations. 

3.2 Steady Flux Limits: We calculated the upper limits to the muon flux from the all-sky survey, 
following the Helene method (Helene, 1983). The steady upper limit at 95% confidence level to the muon flux 
for each bin was computed from: 

rv(95%) 



J^(95%) < 



■cm 2 s 1 



(2) 



Here n M (95%) represents the muon number we would observe at 95% C.L. in each bin, if N j, s (i) and N exp (i) 
are the number of events observed and expected in that bin, respectively. A e j f is the average effective area for 
every bin computed by averaging the projected area seen by each muon 
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where a,i and Si are the muon arrival directions. The A e ff calculation takes into account the geometrical and 
tracking reconstruction efficiencies. T exp is the exposure time, computed bin-by-bin and the factor K takes 
into account the scattering of some muons out of the bin. We found K=0.78 for the selected bin dimension 
(1.5° half-angle). 

We found, for the majority of the bins, J^(95%) < 5x 10" 



1 13 cm 2 s 1 . 



3.3 Search for Modulated Signals: Searches were also made for modulated signals coming from 
those sources that in the past showed variability (Cyg X-3 and Her X-l). We found no evidence for excesses 
in any of the phase bins into which the characteristic period was divided. The 95% C.L. upper limits to the 
modulated muon flux from the directions of Cyg X-3 and Her X-l are, respectively, 1.8 x 10 _13 cm~ 2 s _1 and 
2.2xl(T 13 cm' 2 s- 1 . 

4 Point-like d.c. Source Study: 

To study muons com- 
ing from specific celes- 
tial objects, we selected 
only events contained in 
a narrow cone (1.0° half- 
angle) around the source 
direction. The selected 
windows are centered on 
the Cyg X-3, Mrk 421 
and Mrk 501 positions. 

Table 2 gives the muon Table 2: Number of events observed and expected, the standard deviation calcu- 
number observed and ex- lated as in (I), the average effective area and exposure time for three selected 
pected from the Monte sources. For these we give also the flux upper limits calculated with the Helene 
Carlo simulation. Since formula (Flux 1) and with the classical method (Flux 2). 
there is no significant ex- 
cess from Cygnus X-3, we calculated its steady flux limit following the Helene formula (Helene, 1983) and 
the classical method (Hikasa et al., Particle Data Group, 1992). 

For Mrk 421 and Mrk 501 we have small steady excesses at the level of 1.2a and 2.0cr respectively. For 
these sources the steady flux limits are given in Table 2. 

5 Search for Bursting Episodes: 

During 1997, surface experiments reported several 7- 
bursts from celestial objects like Mrk 421 and Mrk 501. 
Therefore we have made also a search for pulsed muon sig- 
nals in a narrow window (1° half-angle), around the position 
of these possible sources of U.H.E. photons. We studied burst- 
ing episodes following two different methods. 

First, we searched for daily excesses of muon flux above 
the background: because the daily observed events are few, 
we plotted the cumulative excesses day by day, over a long ob- 
servation period (Clay, Dawson, & Meyhandan, 1994). Fig.l 
shows the progressive accumulation of our excesses above the 
computed background for Mrk 501. Periods are present in 
which excesses accumulated rapidly but after that, the cumu- 
lation rate reduced. 

We observe the cumulated signal is every day above the 
background: fluctuations are contained in 2a for Mrk 421 and 
for Mrk 501. 

For the same sources, we have also searched for statisti- 
cally significant daily excesses of muons, following a second 
method (Padilla et al., 1998). From our complete data set we 
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Figure 1 : Cumulative muon excesses from the 
direction of Mrk 501 (1° half-angle). The se- 
lected time period is from January 1992 until 
January 1999 



computed, day by day, the quantity —logio P, where: 
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(4) 
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7 Jan 93 
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1.88 
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2.6xl0~ 4 
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8 


1.72 


3.16 


6.9xl0~ 4 


5 Dec 98 


8 


1.52 


3.48 


3.3 xlO~ 4 



N i, s is the number of observed muons in a day and N bck is the number of background events in the same 
day, computed by Monte Carlo simulation, a is the ratio of the on-source time to the off-source time (0.04 
for the present analysis) and P represents the probability of observing a burst at least as large as N b s due 
to background fluctuations. We assume that the background has a Poissonian distribution about AT bcfc . If no 
bursts are observed, we expect that the cumulative frequency distribution of P is a power law of index -1. 
In Table 3 we report the statistical parameters for Mrk421 to be compared with the observations from surface 
experiments (Shubnell et al., 1996). For comparison, we created a sky survey by computing the quantity -logio 
P for every bin of the sky map: this distribution has slope -1.1 and the superposed line fits the data well. 

6 Conclusions: 

Since February 1989, the MACRO detector has 
collected a sample of about 38.5 million muons. Us- 
ing this sample, we searched for muon excesses above 
background from every sky direction and from single 
cosmic ray sources. No significant excesses have been 
found from the all-sky survey and we computed the 
upper steady limit to the muon flux. At 95% confi- 
dence level, we found this limit at 5 x 10 -13 cm~ 2 s _1 . 
The search for steady emission from different sources 

does not indicate for Mrk 421 and Mrk 501 a muon Table 3: Statistical parameters for possible excesses 
excess above the background, as the statistical signifi- from Mrk421. N exp = a x N bck and P represents 
cance of the excess is 1.2cr and 2a respectively. the probability per day to observe a burst at least 

as large as N obs , due to background fluctuations (1° 
Referefa^eW 1 ,g^' The total number of observations was 

3660. 
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